Butterfly hysteresis curve is a signature of adiabatic Landau-Zener transition 
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We stress that the so-called butterfly hysteresis curves observed in dynamical magnetization mea- 
surements on systems of low-spin magnetic molecules such as V15 and Ve are a signature of adiabatic 
Landau-Zener transitions rather than that of a phonon bottleneck. We investigate the magnetization 
dynamics analytically with the help of a simple relaxation theory in the basis of the adabatic energy 
levels of the spin 1/2, to a qualitative accordance with experimental observations. In particular, 
reversible behavior is found near zero fleld, the corresponding susceptibility being bounded by the 
equilibrium and adiabatic susceptibilities from below and above, respectively. 

PACS numbers: 75.50.Xx,76.60.Es, 75.45.-|-j 



Magnetic hysteresis curves in crystals of molecular 
magnets with an effective spin S* = 1/2 such as V15 
(Refs. iMl and Vg (Ref. 0) have under some conditions 
the so-called butterfly form that is conventionally con- 
sidered as a signature of the phonon bottleneck^ Similar 
phenomenon has been observed on ferric wheels NaFeg, 
where the ground state changes from 5 = to S" = 1 at 
some magnetic field tSi 

We shall demonstrate that butterfly hysteresis curves 
can qualitatively be explained already by a simple re- 
laxational model taking into account adiabatic Landau- 
Zener transitions at avoided level crossing. Using the 
additional equation for nonequilibrium phonons4 does 
not change the butterfly hysteresis curves qualitatively. 
While bad thermal contact between the crystal and the 
holder can apparently lead to the bottleneck and thus 
to smaller experimentally observed effective relaxation 
rates for the spin, these rates could be reduced for other 
reasons, too. In this case butterfly hysteresis curves arise 
without the bottleneck. Although this fact is known (see, 
e.g., Ref. 0), an appropriate quantitative theory is still 
lacking. 

In Vi5 the zero-field splitting A between the two low- 
lying spin levels is about A/fc^ ~50^80 mK (Refs. 1, 
0), and for the experimental sweep rate dB^/dt = 0.1 
T/s the Landau-Zener parameter e = 7rA^/(2?iu) with 
V = 2Sg^gdBz/dt is about 10^. This means that, in 
the absence of relaxation, the system follows the lowest 
adiabatic energy level E- (see Fig.^, the probability to 
remain in the original spin-down state P = being 
negligibly small. 

In the opposite limit of very fast relaxation between 
the adiabatic levels and the system relaxes up 
the energy, E^ E+ , before crossing the resonance and 
down the energy, Ej- — > E- , after crossing the resonance, 
the level populations n± satisfying the equilibrium con- 
dition 
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FIG. 1: Adiabatic Landau-Zener effect with thermal transi- 
tions. The direction of relaxation changes at the energy bias 
Wha that corresponds to the butterfly crossing of the magne- 
tization curve, mz(t) = mi' 



2SgfiB {dBz/dt)t. 

Changing to the spin- up/down basis allows to calculate 
the reduced magnetization 



W 



(n_ - n+) , 



(3) 



where n_ = 1 and rij^ = far before crossing the reso- 
nance. In the absence of relaxation (adiabatic case) one 
has n_ = 1 and rt+ = at any moment of time, while 
at equilibrium n± satify Eq. ^ and the magnetization 
is given by 
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where 



hivo = E+-E^ = VVF2 + A2 



(1) 



(2) 



and W — 2SgiJ.gBz is the energy bias that in the field- 
sweep experiments is a linear function of time: W = vt = 



In the case fc^T 3> A the equilibrium magnetization sim- 
plifies to m'f'^^ = tanh [W/{2kBT)] , independently of A. 
Note that, whatever the relaxation processes, one has 
TOz = at resonance, W = 0. 

In the intermediate-relaxation regime the time- 
dependent values of n± can be obtained from the master 
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FIG. 2: The population difference n-(t) — n+{t) obtained by 
numerical integration in Eq. Q vs the energy bias W for the 
linear sweep from t = — oo to oo in the limit A/fcsT for 
different values of the reduced sweep rate A of Eq. 1151 . 



equation that can be written in the form 



tanh 



hujQ 

2kBT 



(5) 



Note that in the adiabatic hmit e 3> 1 quantum- 
mechanical transitions between the levels E± can be ne- 
glected, and there is no need to use the full density- 
matrix equation. Using the adiabatic basis in Eq. (O 
allows to describe the whole range of the field sweep W, 
in contrast to the master equation in the basis of the 
eigenstates of Sz of Ref.i^that is applicable for ^ A 
only. 

Whereas the equilibrium solution of Eq. (0) for n_ — n+ 
that is attained in the fast-relaxation limit is even in W 
with a minimum at = 0, the general nonequilibrium 
solution for n_ — n+ is lagging behind n^^'^^ — n^'^^ (see 
Fig. 121) . For crossing the resonance in the positive direc- 
tion, it has a minimum at some Wbc > 0. For W < Wbc 



> n 



(eq) 



,(cq) 



one has n_ 

one has n_ — n+ < n'°"'^ — v}^"^' . This overshoot leads 



, whereas for W > Wbc 
,(oq) 



to the butterfly hysteresis curves for mz{t) defined by 
Eq. ©: At W = Wbc the dynamic magnetization mz(t) 
crosses the equilibrium magnetization curve, 



m-,{tbc) 



(4c), 



(6) 



that is the definition of the butterfly crossing. 

The relaxation rate F in Eq. |(SJ is mainly due to the di- 
rect phonon processes and, in the absence of the phonon 
bottleneck, it has the form 



F = Focoth 



/ hug 



(7) 



Here Fq depends on the details of the spin-phonon 
coupling (for V15 and Ve one of the candidates is 
Dzialoshinskii-Moriya interaction^) . We are not going to 
discuss the details of the spin-phonon interactions here as 



the butterfly curve is quite a general phenomenon. One 
has only to take into account that F depends on time via 
the transition frequency ujq- We set 



Fn = F 



00 



(8) 



where in most cases a = 3. For a better fit to experi- 
ments, one can include an tJo-indcpendent relaxation rate 
in Eq. O , as was done in Ref. ,^ The solution of Eq. Q 
has the form 



n_(t)— n-|_(t) 



where 



-,-rt 



rfi'F(t')er(*')*'tanh 



2kBT 



T(t) ^-J^ dt'T{t'). 



(9) 



(10) 



it is easy to check that limt^±oo [n-{t) — n+i^t)] = 1. The 
strongest deviation of n-{t) — n+{t) from 1 originates 
from such t for which tanh[hLUQ{t') / {2kBT)] in Eq. © 
is small enough compared to one. It is the range where 
thermal transitions become significant. 

In the case A ~ kaT that was realized in experiments 
on Vi5 in Refs. (A/fcs ~ 0.05 K, T = 0.1 0.2 
K) the thermal transitions are important in the vicin- 
ity of the level crossing: W — vt ^ A. As the sys- 
tem spends the time tth ^ A/v kBT/v in this re- 
gion, the case Fqo ^ v/A ^ v/{kBT) corresponds to the 
adiabatic limit, n_ = 1 and n+ = 0, whereas the case 
Too ^ v/A v/(kBT) corresponds to the equilibrium 
situation, n^[t) — n^{t) = n'f^\t) — n^^'^\t). In the in- 
termediate case Fqo ^ v/A ~ v/{kBT) the dynamical 
relaxation effect is mostly pronounced, and the hystere- 
sis curve has a butterfly shape with the crossing with 
the equilibrium magnetization curve at Wbc ~ A ~ fcsT. 
Our numerical results obtained from Eq. and shown 
in Fig. I^l are in a qualitative agreement with experi- 
mental and theoretical results of Ref. [J, except for the 
slowest-sweep curve dHz/dt — 0.0044 T/s. 

The low-temperature case A 3> /cgT is trivial, since 
here there are no thermal transitions and the adiabatic 
solution rt_ = 1 and n+ = is valid for all times, while 
the magnetization follows from Eq. 

In the high-temperature range A <g; /cbT (in exper- 
iments on Ve in Ref. the estimated energy gap is 
A/kB — 0.4 K, whereas the maximal temperature was 
T — 4.2 K) thermal transitions take place in a much 
broader region \W\ ^ Wo — kBT than the quantum- 
mechanical crossing \W\ ~ A. The system spends the 
time 



Uh = kBT/v 



(11) 



in the thermal-transition range. Neglecting the narrow 
region |iy| ~ A and approximating Tiloq = \W\ = v\t\ 
one obtains from Eqs. H10() and (O 



fit)t^'^^(^)\c 



v\tl 
knT 



sign(i), (12) 
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FIG. 3: Magnetization hysteresis curves rriz (f) calculated 
from Eqs. @ and for the parameters of Ref. U 
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in the nearly adiabatic limit \ one obtains (for the 
increasing field) 



and 



n^{tbc) - n+{tbc) = tanh 



Wbc 

2kBT 



= 1 



(16) 



where 



dz 



'i>a(l) 



coth ( - 
2 



- 1 



(18) 



/i = 1.60, /2 = 4.62, h = 18.6, A = 94.4, /g = 576. One 
can see that, in fact, the apparent applicability condition 
of Eqs. and (^Tj) is la/ A < 1 that results in rather 
large values of A for a — i and higher. This is an es- 
sential correction to the a priori estimation of Eq. H14|l . 
With increasing a the validity range of Eqs. p6|l and 
(|17|l becomes narrower than Ia/A<^ 1 since the former 
is the first term of the asymptotic expansion in powers 
of 1/ A that becomes progressively bad with increasing a. 
The physical origin of this difficulty is the following. For 
large powers a the situation becomes inhomogeneous in 
the whole domain of 14^: It is nearly adiabatic for smaller 
W and nearly equilibrium for larger W. 

For A ^ fcsT in the nearly equilibrium limit A <C 1 
one obtains 



Wbc^kBTiq^A) 



l/a 



and 



n^itbc) - n+{tbc) = - (qaA) 



l/a 



(19) 



(20) 



FIG. 4: Butterfly-crossing bias Wtc vs the reduced relaxation 
time A for a — 1..5. The nearly-equilibrium asymptotes of 
Eq. p9|l on the left side and the nearly-adiabatic asymptote 
of Eq. 11611 on the right side are shown by the dashed lines. 



where 



$Q (z) = / d-uM^coth 



(13) 



The most pronounced dynamical relaxation effect leading 
to the butterfly hysteresis curve can be expected to take 
place for 

^{tth)tth = — 1, 



A 



where we have introduced 

?! 

A = 



*o(l)' 



(14) 



(15) 



as a reduced relaxation time or reduced sweep rate and 
$i(l) = 2.06, $2(1) = 1-04, $3(1) = 0.670, $4(1) = 
0.527, $5 (1) = 0.424. 

One can find the crossing point Wbc of the butterfiy 

hysteresis curve ruz (t) — mi°'^^ (t) from Eq. by requir- 



ing n^{t) - n+{t) 



,(oq) 



,(cq) 



{t). For A < ksT 



where = (a/2) <i>Q (l)pa and pa is the solution of the 
transcedental equation 



l/a 



T 1 



dppi/"eP-P°. (21) 



One obtains pi = 0.693, p2 = 0.535, pi = 0.482, 
Pi = 0.455, P5 = 0.438, and qi = 0.712, 92 = 0.557, 
93 = 0.506, 54 = 0.480, 94 = 0.464. Note that the gen- 
eral applicability condition of our approximation replac- 
ing fiLOQ ^ \W\ is Wbc ^ A. This yields the applicability 
condition for Eqs. and (pHl 



knT 



< A < 1. 



(22) 



The full dependence of Wbc on the reduced relaxation 
time A of Eq. I|15|l can be found numerically from Eq. 
©. The results for Wbc/iksT) for different a in the 
limit A ^ fcsT are shown in Fig. 0] Note that at the 
butterfly crossing the difference n_ {t) — n+ {t) attains a 
minimum, according to Eq. |SJl ■ The magnetization value 
at the butterfly crossing rri'z'^\tbc) can be obtained from 
Eq. Q and it is given by 



mz{tb 
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^jwlV^ 



ztanh 
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(23) 
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FIG. 5: Susceptibility at the level crossing vs A for a = 3 and 
different values of A/{kBT). The asymptotes given by Eqs. 
I|27|l and Il28|l are shown by the dashed lines. 



Let us now consider the reduced magnetic suscepti- 
bility X = A{dmz/dW) at the level crossing, W = 0. 
Evidently x depends on the sweep rate and is bounded 
by its equilibrium and adiabatic values: 



Xcq = tanh 



A 



2kBT 



< X < 1 = X, 



ad ■ 



With the help of Eq. one obtains the result 
X = n_(0)-n+(0). 



(24) 



(25) 



It is easy to prove that 7i_(0) — n^{0) is the same for the 
increasing field and decreasing field. Then, Eq. (|25|l im- 
plies that the zero-field susceptibility x does not depend 
of the direction of the field. Accordingly the magneti- 
zation, although different from the equilibrium one, is 



reversible near zero field, in full agreement with the ex- 
perimental finding.^ Note that the theoretical aproach of 
Ref. is not applicable in the vicinity of the Landau- 
Zener crossing. Using Eqs. Q and (|25|l one can express 
X as 



° dtrit)e^('^Hanh(^^ 
-oo V ^ksT 



(26) 



that is plotted vs A/Ia for different values of A/ (kBT) in 
Fig. El In the nearly adiabatic limit 1 the asymptote 
of Eq. ||2nil is 



X=l 



2r(Q + 1) 

$a(l)A 



(27) 



In the nearly equilibrium limit A <C 1, the most interest- 
ing is the asymptote corresponding to A/{kBT) — )■ 0: 



X = (a$„ (1) A)^/" 2-(i+i/")r (1 + 1/a) . 



(28) 



Summarizing, we gave a detailed analytical consider- 
ation of the adiabatic Landau-Zener effect with relax- 
ation and we have shown that this is a minimal model 
to describe experimentally observed butterfly hysteresis 
curves. More complicated relaxation models^ includ- 
ing an additional kinetic equation for nonequilibrium 
phonons in the case of the phonon bottleneck seem to be 
nonessential to describe the butterfly hysteresis curves. 
One can interpret the latter in the framework of the sim- 
plest relaxation theory described above by choosing an 
appropriate relaxation rate F that can be effectively re- 
duced because of a poor thermal contact of the crystal 
with the holder. If the details of the phonon dynamics 
related to the phonon bottleneck are the subject of in- 
vestigation, a special care should be taken to single out 
their effect on the magnetization hysteresis curves. 
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